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1.  INTRODUCTION 


The  prediction  of  the  circular  dichroism  of  the  fundamental  vibrational  transitions  of  a  chiral  molecule 
at  the  harmonic  level  of  approximation  requires  the  implementation  of  (Stephens  1985,  1987) 


Rj  =  Im[<0|pei|l>i*<l|Pmag|0>i]  , 
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R|  is  the  rotational  strength  of  the  fundamental  excitation  of  the  ith  normal  mode  of  energy  limj.  The 
corresponding  normal  coordinates,  Qj,  are  related  to  Cartesian  nuclear  displacement  coordinates 
(X  =  nucleus,  a  =  x,y,z)  via 


^Tuol  52  ^Xa,iQi  • 
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The  electric  and  magnetic  dipole  transition  moments,  whose  scalar  products  determine  rotational  strengths, 
are  in  turn  a  function,  respective,  of  the  tensors,  and  M^.  The  tensors  P^,  which  determine 

electric  dipole  transition  moments,  are  the  well-known  atomic  polar  tensors  (APTs)  (Stephens  et  aL  1990). 

X  , 

The  tensors  which  determine  magnetic  dipole  transition  moments,  are  the  atomic  axial  tensors 
(AATs),  more  recently  introduced  into  the  literature  (Stq)hens  1985,  1987;  Stephens  et  aL  1990). 

Since  the  development  of  equation  (Ic)  for  vibrational  magnetic  dipole  transition  momoits,  a 
substantial  number  of  predictions  of  vibrational  circular  dichroism  (VCD)  spectra  utilizing  equation  (1) 
and  ab  initio  SCF  or  post-SCF  c^culations  of  APTs  and  AATs  have  been  published.  For  some  molecules. 
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predictions  have  been  compared  to  experimental  data  (see,  for  example,  Lowe,  Stephens,  and  Segal  1986; 
Jalkanen  et  al.  1987,  1988,  1989,  1990;  Kawiecki  et  al.  1988,  1991;  Bursi,  Devlin,  and  Stephens  1990; 
Lowe  and  Alper  1988;  Dothe,  Lowe,  and  Alper  1988;  Amos,  Handy,  and  Palmieri  1990;  Malon  et 
al.  1992;  Fablioni  et  al.  1993).  Unfortunately,  these  predictions  have  used  vibrational  force  fields  of 
widely  differing  provenance,  reliability,  and  accuracy.  The  accuracy  of  predicted  VCD  intensities  has 
therefore  been  limited  not  only  by  eirors  in  the  APTs  and  AATs  but  also  by  errors  in  vibrational 
fiequencies  and  normal  coordinates.  Unambiguous  attribution  of  the  origins  of  differences  between 
predicted  and  observed  VCD  spectra  has  thus  been  impossible. 

Accurate  harmonic  vibrational  force  fields,  firequencies,  and  normal  coordinates  require  ab  initio 
calculatitMis  at  the  post-SCF  level.  Of  the  various  approaches  to  the  inclusion  of  correlation  in  the 
calculatiwi  of  harmonic  force  fields,  MpUer-Plesset  second-order  perturbation  theory  (MP2)  has  been  the 
most  widely  used.  When  large  basis  sets  are  used,  MP2  force  fields  typically  yield  harmonic  frequencies 
in  error  by  only  l%-3%  (Simandiras  et  al.  1988).  Anharmonicity  corrections  to  vibrational  frequencies 
are  then  comparable  to  or  larger  than  errors  in  harmonic  frequencies.  Calculations  of  vibrational  spectra 
using  such  force  fields  should  not  be  significantly  limited  in  accuracy  by  residual  errors  in  harmonic 
finequencies.  Of  course,  large  basis  set  MP2  calculations  of  harmonic  force  fields  are  computationally 
demanding.  However,  recent  advantages  in  the  computation  of  MP2  second  derivatives  via  analytical 
derivative  methods  (Stanton,  Gauss,  and  Bartlett  1992;  Gauss,  Stanton,  and  Bartlett  1992;  Tracks  et 
aL  unpublished)*  have  greatly  enlarged  the  range  of  molecules  for  which  such  calculations  are  feasible. 
As  a  result,  large  basis  set  MP2  calculations  for  chiral  molecules  whose  VCD  spectra  have  been  measured 
experimentally  are  now  practicable.  We  here  report  calculations  of  VCD  spectra  which  take  advantage 
of  this  development.  The  following  molecules  have  been  studied: 

D 

^ 

2  0  D  2  0  i 

(*= 

MP2  force  fields  have  been  calculated  using  a  [5s4p2d/3s2p]  (TZ/2P  [Stephens  et  al.  1990])  basis  set. 
APTs  and  AATs  have  also  been  calculated  using  this  basis  set.  The  results  are  compared  to  existing  VCD 
spectra  and  to  rotational  strengflis  obtained  thence.  For  1-3  the  sensitivity  of  the  results  to  further 


*  See  GAUSSIAN  92  user’s  manual. 
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enlargement  of  the  basis  set  to  the  [8s6p3d/6s3p]  (VD/3P  [Stephens  et  al.  1990])  level  has  also  been 
examined. 

The  VD/3P  calculations  for  1  were  reported  in  a  prior  communication  (Stephens  et  aL  1993). 

2.  METHODS 

Ab  initio  calculations  have  been  carried  out  using  Cray-YMP  versions  of  CADPAC  4.0  and  5.0  and 
Cray-YMP  and  Cray-2  versions  of  GAUSSIAN  92.  MP2  Hessians  and  APTs  were  calculated  at  the  TZV2P 
basis  set  level  for  oxirane  (5)  and  cyclopropane  (6)  using  ’conventional’  methods  and  CADPAC  4.0  and 
for  propylene  oxide  (4)  using  semidirect  methods  and  GAUSSIAN  92;  the  latter  were  also  used  for 
calculations  on  5  and  6  at  the  VD/3P  basis  set  level.  SCF  APTs  were  calculated  for  each  molecule  at  aU 
basis  set  levels  using  CADPAC  4.0  or  5.0.  AATs  were  ctdculaied  about  an  origin  O  using  the  distributed 
origin  (DO)  gauge  (Stephens  1987;  Stephens  et  al.  1990)  when 

(<r  •  (i^r  *  ^7  E  •  0) 

y5 

where  is  liie  equilibrium  position  of  nucleus  X,  and  (I„p)  is  the  electronic  AAT  of  nucleus  X 
obtained  with  the  coordinate  origin  at  R®.  The  ’local’  AATs,  (I^)^,  were  calculated  at  the  SCF  level 
using  CADPAC  4.0  or  5.0  via  the  relation  (Stephens  et  al.  1990) 

E  •  (“O 

y5 

where  (I^)"  is  the  electronic  AAT  obtained  with  the  origin  at  0  and  E^(jc)  is  the  electronic  APT 
calculated  using  the  momentum  representation  (Stephens  et  al.  1990).  SCF  and  ’semi-MP2’  DO  gauge 
AATs  were  then  obtained  by  combining  SCF  local  AATs  and  either  SCF  or  MP2  APTs  respectively  using 
equation  (3).  Vibrational  frequencies  and  normal  coordinates  were  obtained  from  the  apprc^riately  mass- 
weighted  Hessian.  Rotational  strengths  were  obtained  from  equation  (1).  Note  that  the  use  of  the  DO 
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gauge  in  the  calculation  of  AATs  guarantees  origin-independent  rotational  strengths  (as  long  as  the  same 
APTs  are  used  in  equation  (3)  as  in  equation  [lb]).  Dipole  strengths,  defined  by 

Di  =  l<OlPelU>il^  .  (5) 

were  simultaneously  calculated  via  equation  (lb).  Calculated  frequencies,  dipole  strengths,  and  rotational 
strengths  have  been  used  to  synthesize  unpolarized  absorption  and  VCD  spectra  using  Lorentzian  band 
shapes  (Kawiecki  et  al.  1988). 

3.  RESULTS  AND  DISCUSSION 

Experimental  solution  unpolarized  absorption  and  VCD  spectra  of  1-4  (Kawiecki  et  al.  1991; 
Freedman  et  al.  1991a,  1991b;  Freedman  et  al.  unpublished)  in  C-H  stretching,  C-D  stretching,  and  mid-lR 
spectral  regions  are  reproduced  in  Figures  1-4.  Mid-IR  spectra  have  not  yet  been  reported  for  3. 

Frequencies,  dipole  strengths,  and  rotational  strengths  calculated  using  TZ/2P  MP2  force  fields  and 
APTs  and  T272P  semi-MP2  AATs  for  1-4  are  given  in  Tables  1-4.  Unpolarized  absorption  and  VCD 
spectra  derived  thence  are  plotted  in  Figures  1-4. 

Comparison  of  calculated  and  experimental  unpolarized  absorption  and  VCD  spectra  leads  to 
unambiguous  assignments  of  the  fundamental  bands  of  1,  2,  and  4  in  the  mid-IR  regions  and  of  1-3  in 
the  C-H  and  C-D  stretching  regions.  The  frequencies,  dipole  strengths,  and  rotational  strengths  reported 
(Kawiecki  et  al.  1991;  Freedman  et  al.  1991a,  1991b;  Freedman  et  al.  unpublished)  for  these  bands  are 
listed  in  Tables  1-4.  Assignment  of  the  C-H  stretching  spectra  of  4  is  less  straightforward.  Only  one  of 
the  fundamentals  can  be  confidently  assigned;  its  experimental  fiequency,  dipole  strength,  and  rotational 
strength  are  given  in  Table  4. 

Analogous  calculations  have  also  been  carried  out  for  1-3  using  the  VD/3P  basis  set  The  results  are 
given  in  Tables  1—3  and  plotted  in  Figures  1—3.  All  VD/3P  results  are  qualitatively  identical  to  those 
obtained  at  the  TZJIV  basis  set  level.  They  require  no  changes  in  the  assignments  of  the  spectra  of  1-3. 
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Figure  1.  Predicted  and  experimental  unoolarized  absorption  and  VCD  spectra  of  1,  Upper  and  middle 
spectra  are  predicted  using  TZ/2P  and  VD/3P  basis  sets  respectively;  half  widths  and  half 
height  of  the  Lorentzian  bands  (V)  are  6.0  cm~^  (mid-IR  reeionl  and  10.0  cm~^ 
(C-H/C-H  stretching  regions').  Lower  spectra  are  the  experimental  spectra  of  Frpiftdman 
et  al.  (1991a)  in  CSU  (mid-lR  remon)  and  (C-H/C-D  stretching  regions'):  note  that 

absorption  at  1.034  and  1.117  ctn~^  is  attributed  to  impurities.  VCD  spectra  are  for  (2S.  3SVlT 
Upper  and  lower  frequency  scales  apply  to  calculated  and  experimental  spectra  respectively. 
Left  and  right  vertical  scales  apply  to  C-H/C-D  stretching  and  mid-IR  spectra  respectively. 
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Figure  2.  Predicted  and  experimental  unoolarized  absorption  and  VCD  spectra  of  2.  Upper  and  middle 
spectra  are  predicted  using  TZ/2P  and  VD/3P  basis  sets  respectively:  v  is  6.0  cm~^  and 
10.0  cm~^  in  die  mid-IR  and  C-H/C-D  stretching  regions  respectively.  Lower  spectra  are  the 
experimental  spectra  of  Freedman  et  al.  (unpublished)  in  CS^  (mid-IR  region)  and  CUQ., 
(C-H/C-D  stretching  regions).  VCD  spectra  are  for  CIS.  2SV2.  Upper  and  lower  fireouencv 
scales  apply  to  calculated  and  experimental  spectra  respectively.  Left  and  right  vertical  scales 
apply  to  C-H/C-D  stretching  and  mid-IR  spectra  respectively. 
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Table  1.  Frequencies,  Dipole  Strengths,  and  Rotational  Strengths  of  Molecule  1“ 


a 

& 

O 

ON  Tt  ^  in  Os  p  CN 

i-ioocsoiori'^Tf^ONVom  + 

-  1  1'-' 

u 

(Li 

U 

CO 

VD/3P 

»-^pr7^ppp<svooNt^'-<m^so 

ONiofnTfasi“<^*^osvdi-^r^»-^Ttc5 

^  ^  1  1  <  cn  1  9^ 

1  1  1  ' 

TZ/2P 

plOTfpppp^,  pi-^i-Hp»-jC'|00 

^*ONOONOsr4ococ^ioor^or4c5 

(S)  ^  ^  1  1  1  ^  CO  1 

QO 

o. 

CO 

Q 

> 

iTt"  ^  ^  ^  ^ 

loodpst^r^^cSinoNvdroodocJo 

1— <  <— <  1  11^1  ^  1  'w'  ^  %»✓ 

W  1  W  O' O  w  O  1  nO  w 

•w  O' 

pTft^cscncnpcO'^tspor^ooo 

r4’^*cs^o6^ovdoo6cood^‘^0 

cu 

s 

H 

inodoNr^c^^*c5inod\dcsodc5oc5 

lll*-'l  <N| 

Q 

Xi 

d. 

X 

a> 

^  p  p  p  cn 

cn  lo  cs  ri  Os  od  rf  \d 

lO  ^  ^  ^  ^ 

o 

13 

CJ 

U 

CO 

VD/3P 

p^ppp^p^^^es^Ococsin 
or^ONin<N^vco^odoc4osod*-I 
m  Tf^  cs  ^Ocsosm^ 

CU 

’^cst^r^^pppp'-HCsjospoom 

odNor^ror4^’inooNodoini-Hc4<--I 

rt  CN  OCSONinTT 

% 

Ph 

cn 

Q 

> 

^  ^  ^  ^-N  z-»s  y— s 

ppp^rpTfr^pt^c^t^p^. 

vovdt^'^’odoodr4vdin^T-Icnr4»“< 

'W' 

P'^pr^Tfpr^ppt-ipoor^ONoo 

^’t^odr^o6c5oN^*'or^i-Jriodfoo 

eo  <Ni  ^ 

CSppppTf»“Jp\Ot^CO>Oi-i^,-i 

in  \d  Tt  00  o  ON  vd  Tf  r4  CO  r4  o 

cN  cs  ^  in^cscnT-H 

»-H  T— 4 

Xi 

3027 

3014 

2252 

2232 

1226 

1109 

1102 

948 

914 

\> 

4 

VD/3P 

^i-Hvom^oovoTfcsPSSSSJQJQ 

cnmcacs^-i^-Hi-H,— i,~iO^ONONoor^vo 

o 

1 

r32i^2o5^SS^PiJ^KSQPp} 

c-  -3  3 

I  S  J 

8  P  ^ 
uu  cj 


to 


9 


Calculated  using  SCF  APTs. 

Calculated  using  MP2  APTs  and  senu-MP2  AATs. 

Calculated  using  SCF  APTs  and  AATs. 

Numbers  in  parentheses  are  obtained  using  TZ/2P  APTs  and  AATs. 


Table  2.  Frequencies,  Dipole  Strengths,  and  Rotational  Strengths  of  Molecule  2* 
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From  Freedman  et  al.  (unpublished);  in  C2CI4  for  C-H  and  C-D  stretching  modes  and  in  CS2  solution  otherwise.  Rotational  strengths  were  not  normalized  to  100%  ee;  ee  for 
(IR,  2R)-2  and  (IS,  2S)-i  were  estimated  to  be  >98%. 


Table  3.  Frequencies,  Dipole  Strengths,  and  Rotational  Strengths  of  Molecule  3' 
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I>efinidons  and  units  as  in  Table  1.  Rotational  strengths  are  for  the  (2S,  3S)-enantiomer. 

From  Freedman  et  al.  (1991b);  in  C2CI4  solution  for  C-H  and  C-D  stretching  modes  and  in  CS2  solution  otherwise.  Rotational  strengths  were  not  normalized  to  100%  ee. 


Table  4.  Frequencies,  Dipole  Strengths,  and  Rotational  Strengths  of  Molecule  ^ 


V 

D 

R 

calc. 

TZ/2P 

exp.** 

calc. 

exp.** 

calc. 

exp.** 

MP2 

TZJ2P 

SCF 

Tzm 

MP2 

T272P 

SCF 

TZ/2P 

3259 

3047 

17.3 

29.8 

29.9 

1.9 

2.2 

3.5 

3193 

— 

25.2 

55.1 

— 

-7.1 

-8.0 

— 

3175 

— 

15.9 

30.7 

— 

16.9 

29.0 

— 

3174 

— 

4.8 

12.3 

— 

-19.3 

— 

3156 

— 

21.5 

34.1 

— 

-1.0 

-1.4 

— 

3087 

— 

14.6 

26.5 

— 

-0.4 

-0.6 

— 

1550 

1500 

10.0 

12.0 

15.1 

-3.5 

-43 

^.0 

1529 

1456 

19.9 

15.4 

16.0 

1.8 

1.4 

— 

1514 

1444 

14.1 

10.2 

14.0 

-1.2 

-1.1 

— 

1458 

1406 

51.4 

46.1 

52.5 

-13.4 

-13.2 

-8.4 

1420 

1369 

11.8 

14.3 

14.7 

-4.4 

-3.3 

1295 

1262'' 

11.4 

14.1 

25.7" 

7.5 

15.6" 

1199 

1165 

2.3 

4.0 

4.2 

-0.7 

-1.8 

-0.4 

1177 

1143 

13.0 

20.5 

13.1 

8.5 

9.6 

21.4 

1165 

1130 

2.5 

3.6 

13.9 

-0.8 

-15.9 

1130 

1103 

21.6 

34.7 

26.5 

9.3 

8.1 

1054 

1022 

31.5 

36.6 

56.1 

-7.1 

-8.3 

-6.7 

977 

951 

41.9 

66.1 

65.9 

28.1 

35.5 

71.9 

916 

893 

5.8 

15.0 

19.4 

-17.1 

-24.7 

-42.8 

852 

829" 

199.2 

312.6 

287.6" 

-3.1 

-2.2 

771 

748" 

28.5 

46.8 

40.9" 

-11.7 

-15.2 

-21.6" 

407 

414 

42.7 

52.3 

3.8 

4.5 

371 

373 

46.4 

59.9 

11.9 

15.4 

219 

200'* 

8.1 

10.6 

-3.5 

-3.7 

^  Definitions  and  units  as  in  Table  1.  Rotational  strengths  are  for  the  (S)  enantiomer. 

From  Kawiecld  et  al.  (1991);  in  CCI4  solution,  except  where  indicated.  Rotational  strengths  were  not 
normalized  to  100%  ee. 

^  CS2  solution  data. 

^  Gas>phase  data. 


12 


Calculated  TZ/2P  vibrational  frequencies  arc  in  excellent  qualitative  agrcement  with  experiment  for 
all  assigned  fundamentals  of  1-4.  Calculated  frequencies  are  greater  than  experimental  frequencies,  with 
the  exception  of  two  low  frequency  modes  of  4.  The  ranges/averages  of  the  (absolute)  differences  arc: 
1)  1.5%-6.8%/4.1%,  2)  and  3)  3.1%-7.2%/5.0%,  and  4)  0.5%-7.0%/3.6%.  For  the  mid-IR  modes  alone 
the  ranges/averages  arc:  1)  1.5%-3.4%/2.6%,  2)  and  3)  3.1%-5.0%/3.9%,  and  4)  0.5%-5.0%/3.2%.  For 
the  C-H  and  C-D  stretching  modes  alone,  die  ranges/averages  arc:  1)  5.3%-6.8%/6.1%,  2)  and 
3)  5.5%-7.2%/6.3%,  and  4)  7.0%-7.0%.  Calculated  VD/3P  MP2  frequencies  for  1-3  are  very  close  to 
the  TZ/2PMP2  frequencies.  The  maximum/average  differences  are:  1)  1.2%A).3%,2)  and  3)  2.2%A).7%. 
The  differences  between  calculated  and  experimental  frequencies  arc  very  little  altered.  For  1-3,  a  very 
large  («50%)  increase  in  the  basis  set  size  above  the  T^2P  level  thus  yields  insignificant  improvement. 

The  differences  between  experimental  and  calculation  frequencies  can  be  attributed  to  residual  errors 
in  the  calculated  harmonic  frequencies  and  to  the  contributions  of  anharmonicity  and  solvent  effects  to  the 
experimental  frequencies.  A  comparison  of  calculated  T2y2P  MP2  and  eiq)erimentally  derived  harmonic 
frequencies  for  several  small  molecules  found  differences  averaging  2.2%  (Simandiras  et  al.  1988).  Gas 
phase  and  solution  frequencies  reported  for  1  (Freedman  et  al.  1991a)  and  2  (Freedman  et  al.  unpublished) 
differ  on  average  by  0.4%  and  0.6%  respectively.  We  therefore  conclude  that  anharmonicity  is  overall 
the  dominant  contributor  to  the  differences  in  calculated  and  experimental  frequencies  for  1-4.  This 
conclusion  is  craisistent  with  known  anharmonicity  contributions  to  experimental  frequencies  for  other 
small  molecules  (Hehrc  et  al.  1986). 

Ihpole  strengths  calculated  at  the  TZy2P  basis  set  level  are  in  excellent  qualitative  agreement  with 
experimental  dipole  strengths  for  all  assigned  fundamentals  of  1-4,  demonstrating  the  substantial  accuracy 
of  both  the  hannonic  normal  coordinates  and  the  ATPs.  In  1-3,  dipole  strengths  are  very  little  changed 
on  expansion  of  the  basis  set  to  the  VD^P  level,  indicating  substantial  basis  set  saturation  in  the  TZ/2P 
APTs. 

Quantitative  comparison  of  calculated  and  experimental  dipole  strengths  is  complicated  by  the  absence 
of  weU-defined  error  bars  on  experimental  dipole  strengths,  which  were  obtained  by  fitting  of  absorption 
spectra  using  analytical  band  shape  functions.  The  dipole  strengths  of  weak  bands  and  of  poorly  resolved 
bands  are  particularly  uncertain.  Undoubtedly,  a  significant  fraction  of  the  difference  between  raif^ii^ited 
and  experimental  dipole  strengths  is  attributable  to  experimental  error.  On  the  other  hand,  errors  in 
calculated  dipole  strengths  arising  from  residual  errors  in  nonnal  coordinates  and  APTs  and  fiom  the 
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absence  of  anhannonidty  and  condensed  phase  efforts  can  all  be  expected  to  contribute  significantly. 
Comparison  of  the  dipole  strengths  of  the  dg-  and  d4-isotopomers  of  1,  calculated  at  the  VD/3P  MP2  level, 
to  accurate  gas  phase  experimental  values  demonstrated  (Stephens  et  al.  1993)  that  1)  with  the  exception 
of  C-H  and  C-D  stretching  modes,  calculated  dipole  strengths  agreed  with  experimental  values  within 
experimental  error  while,  2)  in  the  case  of  C-H  and  C-D  stretching  modes,  agreement  was  significantly 
worse.  Errors  arising  from  errors  in  the  harmonic  force  field  and  APTs  and  from  the  absence  of 
anhannonicity  thus  appear  to  be  small  except  in  hydrogenic  stretching  modes,  where  anhaimonicity  is 
more  significant  Given  the  small  differences  in  dipole  strengths  of  1-3  at  the  TZJ2V  and  VD/3P  levels, 
it  can  then  be  inferred  that  for  1-4  at  both  TZ/2P  and  VD/3P  levels:  1)  in  the  mid-IR  region  the 
dominant  contributions  to  the  differences  in  calculated  and  experimental  dipole  strengths  are  condensed 
phase  effects  and  experimental  errors  and  2)  anhannonicity  is  also  likely  to  contribute  significantly  in  the 
C-H  and  C-D  stretching  regions. 

Rotational  strengths  calculated  at  the  TZ/2P  basis  set  level  are  in  excellent  qualitative  agreement  with 
experimental  rotational  strengths  for  all  assigned  fundamentals  of  1-4.  In  particular,  all  signs  are  correctly 
predicted.  The  substantial  accuracy  of  the  MP2  force  field  and  APTs  is  further  confirmed.  In  addition, 
the  substantial  accuracy  of  the  semi-MP2  AATs  is  demonstrated.  In  1-3,  rotational  strengths  are  very  little 
changed  on  expansion  of  the  basis  set  to  the  VD/3P  level,  with  the  exception  of  the  C-H  and  C-D 
stretching  modes.  For  the  mid-IR  modes,  basis  set  saturation  appears  to  be  comparable  to  that  of  the 
dipole  strengths;  for  die  C-H  and  C-D  stretching  modes  it  is  somewhat  less.  The  difference  must  be 
attributed  to  a  lesser  degree  of  basis  set  saturation  of  the  AATs  than  of  the  APTs.  This  conclusion  is 
confirmed  by  calculations  of  dipole  and  rotational  strengths  using  VD/3P  force  fields  and  TZ/2P  APTs 
and  AATs  (Tables  1—3). 

As  are  dipole  strengths,  experimental  rotational  strengths  are  obtained  by  fitting  of  VCD  spectra  using 
analytical  band  shape  functions.  However,  the  quantitative  accuracy  of  experimental  rotational  strragths 
is  substantially  less  well-defined  than  that  of  dipole  strengths  since  experimental  VCD  spectra  arc  less 
accurate  dian  unpolarized  absorption  spectra  (Nafie  1988;  Keideiling  1990).  Resolution  is  lower.  Signal- 
to-noise  ratio  is  much  worse.  Absolute  calibration  is  less  straightforward  and  reliable.  Artifactual  signals 
are  difGcult  to  eliminate.  It  is  consequently  very  difficult  to  quantify  the  contributions  of  experimental 
error  to  the  differences  between  calculated  and  experimental  rotational  strengths. 
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Theoretical  errors  are  also  hard  to  quantitate.  Errors  in  the  harmonic  force  field  and  AFTs  are 
undoubtedly  insignificant  relative  to  errors  in  the  AATs.  The  predominant  deficiency  in  the  AATs  is  the 
incomplete  inclusion  of  correlation  as  a  result  of  the  calculation  of  the  local  AATs  at  only  the  SCF  level. 
In  the  case  of  the  T^P  calculations,  as  discussed  previously,  some  basis  set  error  is  also  present  An 
estimate  of  the  overall  sensitivity  of  rotational  strengths  to  the  inclusion  of  correlation  in  the  AATs  can 
be  obtained  by  eliminating  correlation  contributions  to  the  APTs.  Rotational  strengths  obtained  using  SCF 
APTs  and  AATs  are  given  in  Tables  1-4.  The  (absolute)  changes  in  rotational  strengths  for  those  bands 
of  1-4  for  which  experimental  rotational  strengths  are  available  average  2.9  x  10“^  esu^  cm^  at  the  TZ/2P 
basis  set  level.  The  differences  between  rotational  strengths  calculated  using  T272P  MP2  APTs  and  semi- 
MP2  AATs  and  experimental  rotational  strengths  average  6.4  x  10“^  esu^  cm^.  It  is  therefore  reasonable 
to  conclude  that  a  significant  firaction  of  residual  differences  between  calculated  and  experimental 
rotational  strengths  can  be  attributed  to  the  incomplete  inclusion  of  correlation  in  the  calculation  of  AATs. 

Our  calculations  do  not  successfully  reproduce  the  C-H  stretching  VCD  of  4.  This  could  be  attributed 
to  residual  errors  in  our  harmonic  calculations  or  to  anharmonicity  (Fermi  resonance).  In  view  of  the 
obvious  complexity  of  the  C-H  stretching  absorption  spectrum,  the  latter  is  almost  certainly  the  principal 
contributor. 

A  number  of  calculations  of  rotational  strengths  for  1-4  have  previously  been  reported  (Lowe, 
Stephens,  and  Segal  1986;  Jalkanen  et  al.  1988, 1989, 1990;  Kawiecki  et  al.  1988, 1991;  Lowe  and  Alper 
1988;  Fablioni  1993;  Freedman  et  al.  1991b;  Stephens,  Jalkanen,  and  Kawiecki  1990;  Dutler  and  Rank 
1989;  Rauk  and  Yang  1992;  Yang  and  Rauk  1992;  Polavarapu  and  Bose  1990, 1991;  Freedman  and  Nafie 
1988).  These  calculations  have  used  a  variety  of  harmonic  force  fields;  ab  initio  SCF  (Fablioni  et 
al.  1993;  Dutler  and  Rauk  1989;  Rauk  and  Yang  1992;  Yang  and  Rauk  1992;  Polavarapu  and  Bose  1990, 
1991),  empirically  scaled  ab  initio  SCF  (scaled  quantum  mechanical  [SQM])  (Lowe,  Steffens,  and  Segal 
1986;  Jalkanen  et  al.  1988, 1989;  Kawiecki  et  al.  1988, 1991;  Lowe  and  Alper  1988;  Stephens,  Jalkanen, 
and  Kawiecki  1990;  Freedman  and  Nafie  1988)  and  empirical  general  valence  (Jalkanen  et  al.  1990; 
Stephens,  Jalkanen,  and  Kawiecki  1990;  Polavart^u  and  Bose  1991).  These  diverse  force  fields  vary 
substantially  in  reliability  and  their  deficiencies  are  reflected  in  rotational  strengths  calculated  thence.  The 
harmonic  force  fields  utilized  in  the  present  calculations  are  substantially  more  reliable  than  any  used 
previously  and  reduce  to  insignificance  the  error  in  rotational  strengths  originating  from  errors  in  harmonic 
force  fields. 
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Previous  osculations  have  Sso  adopted  a  considerable  diversity  of  approaches  to  the  calculation  of 
APTs  and  AATs.  APTs  have  been  calculated  ab  initio  at  the  SCF  level  using  both  derivative  methods 
(numerical  [Lowe,  Stephens,  and  SegS  1986;  Lowe  and  Alper  1988]  and  anSytical  [Jalkanen  et  S.  1988, 
1990;  Kawiecki  et  S.  1988,  1991;  Stejiiens,  Jalkanen,  and  Kawiecki  1990;  Duller  and  Rank  1989;  Rauk 
and  Yang  1992;  Yang  and  Rauk  1992;  Polavarapu  and  Bose  1990, 1991])  and  sum-over-states  methods 
(Jalkanen  et  S.  1989;  Fablioni  et  S.  1993).  APTs  have  also  been  calculated  less  accurately  via  fixed 
partis  diaige  (FPC)  (Kawiecki  et  S.  1991;  Stephens,  Jalkanen,  and  Kawiecki  1990),  and  ’floating  basis 
set  vibronic  coupling’  approaches  (Freedman  and  Nafie  1988).  AATs  have  been  cSculated  ab  initio  at 
the  SCF  level  using  both  derivative  methods  (numericS  [Lowe,  Stephens,  and  SegS  1986;  Lowe  and 
Alper  1988]  and  anSyticS  [Jalkanen  et  S.  1988, 1990;  JCawiecki  et  S.  1988, 1991;  Stephens,  Jalkanen, 
and  Kawiecki  1990])  and  sum-over-states  methods  (Jalkanen  et  S.  1989;  Fablioni  et  S.  1993),  and  using 
both  common  origin  (Lowe,  Stephens,  and  SegS  1986;  Jalkanen  et  S.  1990;  Kawiecki  et  S.  1991;  Lowe 
and  Alper  1988;  Fablioni  et  S.  1993)  and  DO  gauges  (Jalkanen  et  S.  1988,  1989,  1990;  Kawiedd  et 
S.  1988,  1991;  Stephens,  Jalkanen,  and  Kawiecki  1990).  In  addition,  AATs  have  been  rairniaM  less 
accurately  via  FPC  (Kawiecki  et  S.  1991;  Stephens,  Jalkanen,  and  Kawiecki  1990),  APT  (Kawiecki  et 
aL  1991;  Stephens,  Jalkanen,  and  Kawiecki  1990),  ’floating  basis  set  vibronic  coupling’  (Freedman  and 
Nafie  1988),  localized  molecular  oibitS  (LMO)  (Polavarapu  and  Bose  1990;  1991),  and  sub-SCF  sum- 
over-stSes  approaches  (Dutler  and  Rauk  1989;  Rauk  and  Yang  1992;  Yang  and  Rauk  1992).  The  MP2 
APTs  and  semi-MP2  DO  gauge  AATs  used  in  the  present  cSculations  are  substantially  more  accurate  than 
any  used  previously  and  considerably  reduce  the  errors  in  rotations  strengths  originating  in  errors  in  APTs 
and  AATs. 

4.  CONCLUSION 

With  the  exception  of  the  C-H  stretching  region  of  4  which  is  strongly  perturbed  by  Fermi  resonance, 
VCD  spectra  of  1-4  cScSated  using  large  basis  set  MP2  force  fields,  together  with  MP2  APTs  and  semi- 
MP2  AATs,  are  in  excellent  agreement  with  experiment  ResiduS  errors  in  frequencies  are  attributable 
predominantly  to  anharmonicity.  ResiduS  errors  in  VCD  intensities  (rotations  strengths)  are  attributable 
predominantly  to  the  incomplete  inclusion  of  correlation  in  the  AATs,  to  anharmonicity  and  to  condensed 
phase  effects.  (In  the  case  of  cSculations  at  the  T^P  level,  incomplete  basis  set  saturation  may  Sso 
contribute  sigmficantly  in  the  C-H  and  C-D  stretching  region.)  Substantially  more  accurate  cScSations 
require:  1)  more  complete  inclusion  of  correlation  in  the  cScSation  of  AATs,  2)  the  inrlnsipn  of 
anharmonicity,  and  3)  the  inclusion  of  solvent  effects.  The  formalism  for  the  cScSation  of  AATs  at  the 
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multiconfiguration  self-consistent-field  (MCSCF)  level  using  either  conventional  or  field-depradent 
(London/gauge  invariant)  atomic  oibitals  has  recently  been  developed  and  implemented  (Bak  et  al.  1993). 
MCSCF  calculations  of  AATs  for  1-4,  should  soon  be  available  and  will  be  of  great  interest  HopeMly, 
calculation  of  AATs  using  other  post-SCF  methodologies  will  also  become  possible  in  the  near  future. 
The  inclusion  of  anharmonicity  in  calculations  of  VCD  spectra  requires  the  extension  of  equations  (lb) 
and  (Ic)  beyond  the  harmonic  level.  While  well-known  for  electric  dipole  transition  moments,  formalisms 
for  magnetic  dipole  transition  moments  including  anharmonicity  have  not  yet  been  developed.  Thus,  while 
great  progress  has  recently  been  made  in  including  anharmonicity  in  ab  initio  predictions  of  vibrational 
frequencies  and  electric  dipole  intensities,  analogous  calculations  of  VCD  spectra  are  not  yet  on  the 
horizon.  The  study  of  solvent  effects  on  vibrational  frequencies  and  intensities  has  a  long  history.  Ihe 
extension  of  current,  state-of-the-art  methods  (macroscopic  and/or  microscopic)  to  VCD  intensities  should 
be  straightforward. 

Our  calctilations  of  the  VCD  spectra  of  1—4  are  substantially  more  accurate  than  previous  work.  For 
the  first  time,  the  harmonic  force  fields  used  yield  frequencies  whose  deviations  from  experimental 
frequencies  are  due  primarily  to  anharmonicity.  The  MP2  APTs  are  substantially  more  accurate  rhan  the 
APTs  used  in  all  earlier  work.  The  combined  use  of  MP2  APTs  and  large  basis  set  SCF  local  AATs  yield 
DO  gauge  AATs  substantially  more  accurate  than  the  AATs  used  in  all  earlier  work.  Comparison  of  the 
present  calculations  and  others  (past  or  future)  which  use  specific  approximations  for  force  fields,  APTs 
and/or  AATs  will  greatly  assist  the  objective  analysis  of  the  relative  accuracies  of  these  approximations. 

The  evaluation  of  the  accuracy  of  our  predicted  VCD  spectra  for  1-4  is  substantially  limited  by  the 
uncertainty  of  the  experimental  VCD  intensities.  This  work  makes  clear  the  importance  of  further 
developments  in  VCD  instrumentation:  reliable  rotational  strengths  of  well-defmed  aceuracy  require 
substantial  improvements  in  resolution,  sensitivity,  absolute  accmracy,  and  artifact  control.  It  is  to  be 
hoped  that  the  greatly  improved  capabilities  of  theory  illustrated  here  will  encourage  such  improvements. 
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NO.  OF 

COPIES  ORGANIZATION 


APPLIED  COMBUSTION  TECHNOLOGY  INC 
ATTN  A  M  VARNEY 
P  O  BOX  607885 
ORLANDO  FL  32860 

APPLIED  MECHANICS  REVIEWS 
ASME 

ATTN  R  E  WHITE  &  A  B  WENZEL 
345  E  47TH  STREET 
NEW  YORK  NY  10017 

TEXTRON  DEFENSE  SYSTEMS 
ATTN  A  PATRICK 
2385  REVERE  BEACH  PARKWAY 
EVERETT  MA  02149-5900 

BATTELLE 

TWSUAC 

505  KING  AVENUE 

COLUMBUS  OH  43201-2693 

COHEN  PROFESSIONAL  SERVICES 
ATTN  N  S  COHEN 
141  CHANNING  STREET 
REDLANDS  CA  92373 

EXXON  RESEARCH  &  ENG  CO 
ATTN  A  DEAN 
ROUTE  22E 
ANNANDALE  NJ  08801 

GENERAL  APPLIED  SCIENCE  LABS  INC 
77  RAYNOR  AVENUE 
RONKONKAMA  NY  11779-6649 

GENERAL  ELECTRIC  ORDNANCE  SYSTEMS 
ATTN  JMANDZY 
100  PLASTICS  AVENUE 
PITTSFIELD  MA  01203 

GENERAL  MOTORS  RSCH  LABS 
PHYSICAL  CHEMISTRY  DEPARTMENT 
ATTN  T  SLOANE 
WARREN  MI  48090-9055 

HERCULES  INC 
ATTN  W  B  WALKUP 
E  A  YOUNT 
P  O  BOX  210 

ROCKET  CENTER  WV  26726 


1  HERCULES  INC 

ATTN  R  V  CARTWRIGHT 
100  HOWARD  BLVD 
KENVIL  NJ  07847 

1  ALLIANT  TECHSYSTEMS  INC 

MARINE  SYSTEMS  GROUP 
ATTN  D  E  BRODEN  MS  MN50  2000 
600  2ND  STREET  NE 
HOPKINS  MN  55343 

1  ALLIANT  TECHSYSTEMS  INC 

ATTN  R  E  TOMPKINS 
MN  11  2720 

600  SECOND  ST  NORTH 
HOPKINS  MN  55343 

1  IBM  CORPORATION 

RESEARCH  DIVISION 
ATTN  A  C  TAM 
5600  COTTLE  ROAD 
SAN  JOSE  CA  95193 

1  nr  RESEARCH  INSTTIUTE 

ATTN  R  F  REMALY 
10  WEST  35TH  STREET 
CHICAGO  IL  60616 

1  LOCKHEED  MISSILES  &  SPACE  CO 

ATTN  GEORGE  LO 
3251  HANOVER  STREET 
DEPT  52  35  B204  2 
PALO  ALTO  CA  94304 

1  OUN  ORDNANCE 

ATTN  V  MCDONALD  LIBRARY 
P  O  BOX  222 
ST  MARKS  FL  32355-0222 

1  PAUL  GOUGH  ASSOCIATES  INC 

ATTN  P  S  GOUGH 
1048  SOUTH  STREET 
PORTSMOUTH  NH  03801-5423 

1  HUGHES  AIRCRAFT  COMPANY 

ATTN  T  E  WARD 
PO  BOX  11337 
TUCSON  AZ  85734-1337 


25 


NO.  OF 
COPIES 

1 

3 


1 

1 

1 

3 


3 

1 

1 


ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


SCIENCE  APPLICATIONS  INC 
ATTN  R  B  EDELMAN 
23146  CUMORAH  CREST 
WOODLAND  HILLS  CA  91364 

SRI  INTERNATIONAL 
ATTN  G  SMITH 
D  CROSLEY 
D GOLDEN 

333  RAYENSWOOD  AVENUE 
MENLO  PARK  CA  94025 

STEVENS  INSTITUTE  OF  TECH 
DAVIDSON  LABORATORY 
ATTN  R  MCALEVY  m 
HOBOKEN  NJ  07030 

SVERDRUP  TECHNOLOGY  INC 
LERC  GROUP 

ATTN  R  J  LOCKE  MS  SVR  2 
2001  AEROSPACE  PARKWAY 
BROOK  PARK  OH  44142 

SVERDRUP  TECHNOLOGY  INC 
ATTN  J  DEUR 

2001  AEROSPACE  PARKWAY 
BROOK  PARK  OH  44142 

THIOKOL  CORPORATION 
ELKTON  DIVISION 
ATTN  R  BIDDLE 
R  WILIER 
TECH  LIB 
P  O  BOX  241 
ELKTON  MD  21921 

THIOKOL  CORPORATION 
WASATCH  DIVISION 
ATTN  S  J  BENNETT 
P  O  BOX  524 
BRIGHAM  CITY  UT  84302 

UNITED  TECHNOLOGIES  RSCH  CENTER 
ATTN  A  C  ECKBRETH 
EAST  HARTFORD  CT  06108 

UNITED  TECHNOLOGIES  CORP 
CHEMICAL  SYSTEMS  DIVISION 
ATTN  R  R  MILLER 
P  O  BOX  49028 
SAN  JOSE  CA  95161-9028 


1  UNIVERSAL  PROPULSION  COMPANY 

ATTN  H  J  MCSPADDEN 
25401  NORTH  CENTRAL  AVENUE 
PHOENIX  AZ  85027-7837 

1  VERTTAY  TECHNOLOGY  INC 

ATTN  E  B  HSHER  ' 

4845  MILLERSPORT  HIGHWAY 
EAST  AMHERST  NY  14051-0305 

1  FREEDMAN  ASSOCIATES 

ATTN  E  FREEDMAN 
2411  DIANA  ROAD 
BALTIMORE  MD  21209-1525 

3  ALLIANT  TECHSYSTEMS 

ATTN  C  CANDLAND 
L  OSGOOD 
R BECKER 
600  SECOND  ST  NE 
HOPKINS  MN  55343 

1  US  ARMY  BENET  LABORATORY 

ATTN  SAM  SOPOK 
SMCAR  CCB  B 
WATERVLIET  NY  12189 
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NO.  OF 

COPIES  ORGANIZATION 


ABERDEEN  PROVING  GROUND 


36  DIR  USARL 

ATTN:  AMSRL-WT-P,  A  HORST 
AMSRL-WT-PC, 

R  AFIFER 
GF  ADAMS 
W  R  ANDERSON 
R  A  BEYER 
S  WBUNTE 
C  F  CHABALOWSKI 
K  P  MCNEBLL-BOONSTOPPEL 
A  COHEN 
RCUMPTON 
R  DANIEL 
D  DEVYNCK 
NFFELL 
B  E  FORCH 
JMHEIMERL 
A  J  KOTLAR 
MRMANAA 
WFMCBRATNEY 
K  L  MCNESBY 
S  VMEDLIN 
MSMILLER 
AWMIZIOLEK 
S  H  MODIANO 
J  B  MORRIS 
J  E  NEWBERRY 
S  A  NEWTON 
R  A  PESCE-RODRIGUEZ 
BMRICE 
R  C  SAUSA 
M  A  SCHROEDER 
J  A  VANDERHOFF 
M  WENSING 
AWHREN 
J  M  WIDDER 
C  WILLIAMSON 
AMSRL-CI-CA,  R  PATEL 


Intentionally  left  blank. 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  tmdertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your  comments/answers 
to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  ARL  Report  Number  ARL-TR-964 _ Date  of  Report  Febmarv  19% _ 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for  which  the  report 

will  be  used.) _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of  ideas,  etc.) 


S.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved,  operating  costs 
avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate  changes  to 
organization,  techiucal  content,  format,  etc.) _ 


Organization 

CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 

City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address  above  and  the 
Old  or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  tape  closed,  and  mail.) 
(DO  NOT  STAPLE) 


